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Abstract 

We consider quantum electrodynamics quantized on the light front in Feyn- 
man gauge and regulated in the ultraviolet by the inclusion of massive, 
negative-metric Pauli-Villars (PV) particles in the Lagrangian. The eigen- 
state of the electron is approximated by a Fock-state expansion truncated 
to include one photon. The Fock-state wave functions are computed from 
the fundamental Hamiltonian eigenvalue problem and used to calculate the 
anomalous magnetic moment, as a point of comparison. Two approaches 
are considered: a sector-dependent parameterization, where the bare pa- 
rameters of the Lagrangian are allowed to depend on the Fock sectors be- 
tween which the particular Hamiltonian term acts, and the standard choice, 
where the bare parameters are the same for all sectors. Both methods are 
shown to require some care with respect to ultraviolet divergences; neither 
method can allow all PV masses to be taken to infinity. In addition, the 
sector-dependent approach suffers from an infrared divergence that requires 
a nonzero photon mass; due to complications associated with this divergence, 
the standard parameterization is to be preferred. We also show that the self- 
energy effects obtained from a two-photon truncation are enough to bring the 
standard-parameterization result for the anomalous moment into agreement 
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with experiment within numerical errors. This continues the development of 
a method for the nonperturbative solution of strongly coupled theories, in 
particular quantum chromo dynamics. 
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1. Introduction 

The nonperturbative solution of quantum field theories has proven to 
be a difficult task. The method that has had the most success to date, 
lattice gauge theory jlj, has attained this success only after a long period of 
development, with a great number of technical innovations along the way. 
What is more, the level of success that can be achieved is inherently limited 
by the lack of direct contact with wave functions for bound-state constituents 
and by the Euclidean formulation. The limitation of Euclidean formulation is 
shared by the method of Dyson-Schwinger equations [2j , which are coupled 
equations for the n-point Euclidean Green's functions where bound states 
appear as poles in the propagators. Solution of the infinite system requires 
truncation and a model for the highest n-point functions. 

In order to calculate wave functions directly in a Minkowskian formula- 
tion, there has been an effort for a number of years to develop a Hamiltonian 
approach in light-cone quantization j^]. Although the ultimate objective is 
to be able to solve for the bound states of QCD, most of the development 
has been in QED and simpler theories, such as Yukawa theory and <f) A the- 
ory. Many two-dimensional theories have been solved; however, success with 
four- dimensional theories has been limited by the need for a consistent regu- 
larization and renormalization scheme and by the large size of the numerical 
calculation. There is also a related light-front lattice Hamiltonian formula- 
tion, known as the transverse lattice method j4|], which we do not consider 
here. 

A regularization scheme that has been useful in doing perturbative light- 
cone calculations is the alternate denominator method of Brodsky, Roskies, 
and Suaya 0] . Because of its success, one would naturally consider extending 
this approach to nonperturbative calculations. Unfortunately, application of 
the method requires explicit identification of light-cone energy denomina- 
tors, which are only implicit in the coupled equations of the nonperturbative 
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mass eigenvalue problem. One could instead use the alternate denominator 
method to construct counterterms, and then incorporate the counterterms in 
the Hamiltonian for the nonperturbative calculation. However, these coun- 
terterms would be limited to a particular order in the coupling while the 
nonperturbative problem sums a partial set of contributions to all orders. 
Instead, one needs a method that generates counterterms to all orders, but 
no more than what is needed. 

A method for regularization that has proven quite useful is Pauli-Villars 
(PV) regularization J6f , which was developed and tested in a series of calcu- 
lations 

a a a 

[lol . Ell . Ilil . [lil . Q. The key idea is to include enough PV 



fields in the Lagrangian to regulate the theory perturbatively and, where 
possible, maintain symmetriesjj The derived light-front Hamiltonian then 
defines the nonperturbative bound-state problem. As for the eigenstate, it 
is approximated by a truncated Fock-state expansion. Then the mass eigen- 
value problem leads to coupled integral equations for the Fock-state wave 
functions. The PV particles appear in the Fock states and, through chosen 
negative metrics, bring about the subtractions necessary to regulate the in- 
tegral equations. A possible formulation for QCD along these lines has been 
given by Paston et al. (l6| . 

Two methods of parameterization are in use. One is the standard choice, 
where the bare parameters are those of the regulated Lagrangian. The other 
is a sector-dependent parameterization [l7j], where the bare parameters are 
allowed to depend on the Fock sector(s) on which the Hamiltonian acts. In 
either case, the parameters are fixed by constraints from observables and from 
symmetry restorations. The standard scheme has been used extensively in 
studies of PV regularization 10, Ilil, ll2l, |13|, |14| . The sector-dependent 



scheme was first systematically applied to QED by Hiller and Brodsky [18 



though they did not consider a sector-dependent vertex mass. More recently, 



it was investigated by Karmanov, Mathiot, and Smirnov [19] . In order to 
better understand how to proceed with the large-scale numerical calculations 
that need to be done, a comparison of these approaches needs to be made. 

Both parameterizations require some care, particularly with respect to 
uncanceled divergences. In the standard method, an uncanceled divergence 



appears in the following way [12[. The results for a generic physical quantity, 



1 The introduction of PV partners to the fields of a theory has recently been used to 
define extensions of the Standard Model that offer a solution to the hierarchy problem (HI . 
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such as the electron's anomalous moment, will be of the form 



„ 2 a i [+ a 25 ,2 ln/ipv H ] f 0, with truncation 

1 + £»! t7 2 + b 2 g 2 In /x PV + ■ ■ • \ finite, without truncation, 

(1) 

where /xpy is a PV mass scale and the contents of the square brackets are 
absent in the case of truncation. When the limit /ipy — > oo is taken, the 
result is either zero or a finite value. As the PV mass fipy is increased, 
the divergence in the denominator is not canceled by a divergence in the 
numerator, if the truncation is made. Therefore, any potentially meaningful 
calculation must be done at finite /ipy, in a range of scales where the errors 
due to truncation and PV inclusion are minimized. 

The sector-dependent method is also limited to finite PV mass scales by 



an uncanceled divergence. As shown in [18|], the bare coupling can become 
imaginary if the PV masses become too large. Also, as we will show below, 
the probabilities for individual Fock sectors are driven outside the range of 
to 1. The authors of [19( do not consider these limitations and calculate 
quantities that do not depend directly on this bare coupling; they then have a 
result at infinite PV masses, but the underlying theory and its wave functions 
are ill defined. 

Instead, the calculation should be done at finite PV masses. Here we do 
this, for the anomalous magnetic moment of the electron computed in the 
one-photon truncation of the dressed-electron state, and compare the results 
of the standard [ljj and sector-dependent 19] approaches. In this particular 
context, the dressed electron can be viewed as a bound state of bare con- 
stituents, electron and photon, in QED. The particular truncation is chosen 
because the eigenvalue problem can be solved analytically and the compari- 
son is not obscured by numerical artifacts. The calculation of the anomalous 



moment is not intended as competitive with perturbation theory [20j; the 
weak coupling of QED makes nonperturbative calculations much less accu- 
rate. However, the calculation does stand as a test of a method that can be 
applied where perturbation theory is useless. 

Convergence with respect to the Fock-state expansion has been checked in 
Yukawa theory jsj. There a full discrete light-cone quantization (DLCQ) [2l[ 
calculation was possible. The DLCQ approach does truncate the Fock-state 
expansion but only to the extent required by the numerical resolution; more 
constituents are allowed as the discrete momentum fraction is reduced. Ex- 
plicit truncations to fewer constituents can be made and the results com- 
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pared. For Yukawa theory, the comparison showed that the Fock-state ex- 
pansion converges quickly. 

The analogous results for a cubic scalar theory |22| do not show such 
rapid convergence. However, this is because the spectrum of the theory is 



unbounded from below [23|], and the lowest eigenstates are dominated by 
Fock states with large numbers of constituents (24| . In other words, for 
cubic scalar theories the highest Fock states make the largest contributions. 
Yukawa theory and QED do not suffer this fate because of the Pauli principle. 
In order to have a simple vacuum and a well-defined Fock-state expansion, 



we use light-cone coordinates [25|, |3J x = (x~ = t — z,x± = (x, y)) and light- 
cone momenta p = (p + = p° +p z ,p±). The light-cone energy p~ = p° — p z is 
conjugate to the light-cone time x + = t + z. Stationary states are obtained 
as eigenstates of the light-cone Hamiltonian |3| V~ , which defines the mass 
eigenvalue problem V~\P) = ^+\P), in a frame where the total transverse 

momentum P± is zero. Further details of the coordinate choice can be found 
in 



14| or |3[. 

To further develop the method in the context of a gauge theory, we con- 
sider QED. To regulate QED, we use one PV electron and two PV photons. 
The second PV photon is needed to restore chiral symmetry in the massless- 



electron limit [141 ] . The couplings and metrics of the PV fields are adjusted 
to accomplish the ultraviolet regularization and the finite chiral-symmetry 
correction. They also provide for cancellation of instantaneous fermion in- 
teractions and allow the constraint equation for the nondynamical fermion 



fields to be solved exactly without use of light-cone gauge [14 



In Sec. 2 of [l2j, there is some discussion of the use of three PV fermions 
to regulate QED. However, this was included only to show that naive PV 
regularizations in light-cone gauge do not necessarily work. A regularization 
in light-cone gauge that does work is given in Sec. 4 of that paper, but 
this requires additional regulators that complicate the theory. The simplest 
regularization, which uses one PV fermion and one PV boson, is found for 
Feynman gauge and discussed in Sec. 3. There the chiral limit is correct, but 
only because the PV electron mass is taken to infinity. If the PV electron 
mass is kept finite, regularization by one PV photon is insufficient for both 
standard and sector-dependent parameterizations. 



Also, the result in 121 for the anomalous moment of the electron differs 



from the Schwinger result 26]. The source of the difference is not the un- 



canceled divergence. Instead, it is contributions from all higher orders in the 
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coupling constant, consistent with the Fock-space truncation to one electron 
and one photon. The effect of the uncanceled divergence, which appears in 
the denominators of Fock-state probabilities, is logarithmic and quite mild, 
and, in any case, can only make the estimate for the anomalous moment 
smaller than the Schwinger result. The estimate obtained is instead larger, 
due to sensitivity to the constituent (bare) electron and photon masses and 
the nonperturbative shift of the constituent electron mass. Here we correct 
this shift by including the self-energy contribution from the one-electron/two- 
photon sector and find very good agreement with exp eriment. 

The sector-dependent approach, as used in fl9f, does no better than 
lowest-order perturbation theory, which is quite strange, since a nonper- 
turbative calculation should include some physics to all orders. In addition, 
the calculation in [l9[ suffers from uncanceled divergences that appear in the 
numerators of Fock-state probabilities, instead of the denominators, and that 
push these probabilities outside the physical range of [0, 1]. 

With this as introduction, we continue in Sec. [2] with a comparison of the 
standard and sector-dependent parameterizations in a one-photon truncation 
of the Fock expansion for the dressed electron. In Sec. [31 the standard- 
parameterization results are extended to include self-energy contributions 
that correspond to intermediate states with two photons. Section H] contains 
a summary of the comparison and some general discussion of applications. 



Appendix A provides specifics about light-front QED in Feynman gauge, 



with further details available in [12| and [14j . Appendix B gives details of 



the solution when self-energy contributions are included. 



2. One-Photon Truncation 

We work in Feynman-gauge QED quantized on the light-front and reg- 
ulated with one PV electron and two PV photons. Details are given in 
Appendix A 

2.1. Fock- State Expansion 

We seek the solution for the dressed-electron eigenstate truncated to in- 
clude only the bare-electron and one-electron/one-photon Fock sectors. The 
Fock-state expansion for the eigenstate with total J z = ±| is then 

I^GQ) = J2 ^L(£)|0> + £ / dkC£(k)bt(P - k)al(k)\0}- (2) 
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The amplitudes Zi and wave functions C£ a define this state. 

The interpretation of the expansion requires projection onto a physical 



subspace. We apply the same approach as was used in Yukawa theory [11 
where a projection onto the physical subspace is accomplished by expressing 
Fock states in terms of positively normed creation operators and null com- 
binations. Here, the positive-norm operators are aJ , ajL, and bl s , and the 
null combinations are = V^ !^ an d b\ = 6j s + b\ s . The b\ particles are 
annihilated by the generalized electromagnetic current ?/>7 M ?/>; thus, b\ cre- 
ates unphysical contributions to be dropped, and, by analogy, we also drop 
contributions created by a*. The projected dressed-fermion state is 

K hys (£)) = E(- 1 ) J ^o ± (aio) (3) 

i 

+ £/*££ £ { -^Ct s (k)bl(P-k)al(k)\0)- 

s/i J i=0 3=0,2 k=j/2 V Sfc 

The normalization condition 

(r P U^)\^ys(p)) = - P) 5 °'° ( 4 ) 

is then applied. 

The projection is critical, not only for the removal of negatively normed 
contributions which could make probabilistic interpretations difficult, but 
also for the regularization of expectation values. In particular, a calculation of 
the expectation value of the number of photons in the dressed-electron state 
will yield infinity for the unprojected Fock expansion, because the spin-flip 
wave function C^-(k) falls off inversely with k±. The projection introduces 
subtractions that remove this divergent behavior. In [19j, such a projection 
is not made and instead the eigenstate is to be interpreted only as part of a 
larger process; the associated Fock-state wave functions have no direct utility. 

2.2. Coupled Equations 

Our Fock-space-truncated eigenstate has to satisfy the mass-squared eigen- 
value problem V~\ip) = j£+\?p). Projections onto each Fock sector yield cou- 
pled equations for the one-photon amplitudes Zi and one-electron/one-photon 
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wave functions C^. These coupled equations are, with y = k + /P + , 

[M 2 -m 2 }z t = f(P + ) 2 dyd 2 k ± J2Vli(-^ J+l ^ (5) 



and 



M 2 - ™* + kl 



M 2 m i +kl 



3,1,1* 

x [Vjt* ± (E ~ £, P)Cf ± (k) + Ufe(P - k, P)Cf T {k)] 

y 
y 



ct ± {k) = v / ^E(- 1 ) i '^ p+ ^ ± (£-^a, 

V 

(6) 

%' 

(7) 

The indices are arranged such that an index of i corresponds to the one- 
electron sector and j to the one-electron/one-photon sector. Therefore, in 
the sector-dependent approach, a mass rrii in a vertex function is assigned the 
bare mass, and rrij is the physical mass. In the standard parameterization, 
all are bare masses. 

The coupled equations can be solved analytically 12j . The wave functions 

E,(-i)^ + ^ ± (p-^p) 



C£ are 



eg® 



t ^±3 

M 2 _ m 1+ k l _ tf+ k l 

t E,(-i)^ + ^ ± (Z-i,Z) 



M 2 



}+k 2 tf+k 2 



1-2/ 



and the amplitudes satisfy 

(M 2 - m 2 ) Zi = 2e 2 ^(-l)' 4 '^ [J + m(m v h ~ 2(m* + 



with 12 



I n (M 2 ) 
J(M 2 ) 



dydk] 

16tt 2 ^ M2 _. <";-'--i _ /'r+/-j //{ 1 - u)" 



dydk\ sr-^ 



i-y y 



m 2 + k\ 



16tt 2 



it!i_iM!/(l-2/) 2 ' 

1-1/ 2/ 



M 2 
8 



(8) 
(9) 

(10) 

(11) 
(12) 



Each integral can be computed analytically. For the sector-dependent pa- 
rameterization, the convention of m ; being a bare mass and rrij a physical 
mass has been maintained, with the extension that rrii> is also a bare mass. 
The integrals I n and J then depend only on the physical mass, and fflOl) is 
equivalent to perturbation theory. For the standard parameterization, I n and 
J depend on the bare mass, and (1 1 01) is nonperturbative. 



With use of the identity 14j J = M 2 I , J can be eliminated from the 



eigenvalue problem. This allows the solution to the eigenvalue problem (flOl) 



to take the simple form [12 



(M±m )(M±m 1 ) ^ = M±m ^ 
Qo± 87r(mi-m )(2J 1 ±MJ )' Zl M± mi Z °' [ } 

with a = el/4iT and z determined by the normalization (J3J), which reduces 
to 

1 = (z - Zl ) 2 (14) 

"* 1,1' i'i j=0,2 

i/2+li/2+l , , k , +k 

E E 



; .. 2 \yrnl + {l-y)nl + kl-mly{l-y)\ 



m v mj - (rrii + mj/)(m; + m;/)(l - y) + mim v {\ - y) 2 + 
[ymf + (1 - y)n\ + k 2 ± - m\y{\ - y)\ 

2.3. Discussion of the Solution 

Since the PV fermion is needed only to regulate the integral J, which 
has been eliminated from the calculation, we can safely take the m\ — > oo 
limit, to simplify the remaining steps. In this limit, we have z\ = 0, ni\Z\ — > 
±(M — m )zo, and 

M(M±m ) 

" 0± = ± 8 7 r(27 1 ±M/ o) - (15) 



Also, the third constraint ( 1A.5[) is automatically satisfied, and the second PV 



photon flavor can be discarded. 

For the standard parameterization, we cannot solve explicitly for mo, 
because Iq and l\ are functions of mo- However, «o± is just ct± and the value 
of m is determined by requiring a>± to be equal to the physical value of a. 



9 



This defines a nonlinear equation for m , due to the dependence of I n on m . 
For small values of the PV masses there may be no such solution; however, 
for reasonable values we do find at least one solution for each branch. The 



plot in Fig. 1 of [14J shows a±/ct as functions of mo- The a_ branch is the 
physical choice, because the no-interaction limit (a_ = 0) corresponds to 
the bare mass becoming equal to the physical electron mass, M = m e . 
This is consistent with the sector-dependent case, where the lower sign is also 
chosen. 

If the PV electron has a sufficiently large mass, the value of mo that yields 
a_ = a is less than m e . In this case, the integrals I n and J contain poles 



for j = I = and are defined by a principal- value prescription [121 ] . For 
terms in the normalization sum fl!4p with % = k = or i' = k! = 0, there are 
simple poles, again defined by a principal-value prescription. For the terms 
where all four of these indices are zero, there is a double pole, defined by the 
prescription [12] 



dx . ^ N = lim — 
(x — a) 1 n^o 2i] 



V I dx^ vldx- f{x) 



x — a — 7] J x — a + rj 



(16) 



In the sector-dependent approach, Ji and Jo are independent of m , and 
the solution for «o in (|T5|) can be written as an explicit expression for m 

m = T M + 87^(2/! ± M/ Q ). (17) 

The lower sign is chosen, so that mo reduces to M = m e when the coupling 



is zero. This expression for mo is equivalent to Eq. (74) of 19], since the 
mi — > oo limits of Io and I± are equivalent to B and A in (lif . 

There is, however, the sector dependence of the coupling. For the cou- 
pling between the bare-electron and one-electron/one-photon sectors, the 
bare coupling is given by Eq. (3.20) of [l8|, which in our notation is writ- 
ten e = e/zo, where z is the amplitude for the bare-electron Fock state 
computed without projection onto the physical subspace. This expression 
arises from the following considerations. In general, the bare coupling would 



be eo = Zie/ Z2iZ 2 fZ 3 ; this includes the truncation effect that splits the 
usual Z2 into a product of different \fZ~2 from each fermion leg [27]] . Since 
no fermion- ant if ermion loop is included, we have Z s = 1. Since only one 
photon is included, there is no vertex correction and Z\ — 1. Also, only the 
fermion leg with no photon spectator will be corrected by \fZ~2 and therefore 

\fZ2iZ2~f = Zq. 
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In the infinite-mx limit, the bare-electron amplitude without projection 
is determined by the normalization 

1 = z 2 + elz 2 Q J 2 , (18) 

obtained as a limit of the unprojected form of ([HP , with 

On replacement of eo by e/zo, we can solve for zq as 



z = ^l-e 2 J 2 (20) 

and also find 



e = e/yi-e 2 J 2 . (21) 

This result matches Eq. (79) of [19|], with J 2 in our Eq. (TT9"j) agreeing with 
Eq. (77) in 19J. However, our interpretation is quite different, in that we do 
not allow eo to become imaginary, which limits J 2 . Since J 2 grows logarith- 
mically with fii, as 

1 ^ + (22) 



87H \ m„ 



this last PV mass cannot be taken to infinity. In [19], this behavior is dis- 
regarded, /ii is taken to infinity, and the eigenstate is re-interpreted in the 
context of some larger process, where eo is not directly referenced. Un- 
fortunately, an imaginary eo makes the underlying theory unphysical, with 
opposite charges repelling each other. There would also be strange results 
for probabilities when the uncanceled divergence is ignored, so that J 2 grows 
without limit. The probability for the one-electron sector is z\ = 1 — e 2 J 2 < 
and for the one-electron/one-photon sector, e 2 J 2 > 1. 

For the standard parameterization, there is also a limit on The pro- 
jected normalization condition ( [14]) can be written as 

1 = zl + e 2 zl J 2 , (23) 

with 

■h = r-j / V dydkl[m 2 - 4m m e (l -y) + m 2 e (l- y) 2 + k]_\ (24) 




2 



[k± + (1 - y)l4 + yml - y{l - y)m 2 e ] 
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Thus the bare amplitude is 



^ = l/Vl + eV 2 , (25) 

which is driven to zero as fix — > oo and causes most expectation values also to 
go to zero. This is again the difficulty of uncanceled divergences 12J. Here the 
denominator contains a term one order higher in a than the truncation allows 
in the numerator. Any divergence in this higher-order term is not canceled 
by a corresponding divergence in the numerator, and the expectation value 
becomes zero. 

To obtain meaningful results for both standard and sector-dependent ap- 
proaches, one must balance minimization of the truncation error against min- 
imization of the errors associated with having finite-mass, negative-metric PV 



particles in the basis [12j. We look for a range of PV masses where expec- 
tation values are approximately constant, and find that this does occur for 
the anomalous moment. Since J 2 and J 2 grow logarithmically with fii, the 
dependence is very mild, and H\ can be taken quite large. 

There is an additional complication of an infrared divergence in the sector- 
dependent case. The pole in the wave function (jS]) moves to an endpoint 
when fi = and the constituent mass m is set equal to the physical mass 
m e . This happens because the denominator of the wave function at zero 
transverse momentum is proportional to m\ — m^/(l — y) — ^/y. At the 
endpoint, a principal- value prescription cannot be used, and the mass of 
the physical photon hq must be nonzero to remove this pole. However, the 
anomalous moment of the electron is particularly sensitive to this mass, and 



even a small value for /i will cause a noticeable shift 18]. From another 
side, the integral J2 is driven to negative values as /j,q goes to zero, as can 
be seen from (1221) . making the probability of the two-particle sector negative 
and therefore unphysical. Thus, the fi — > limit is needed for an accurate 
result, but the limit cannot be taken. The standard parameterization does 
not suffer from this infrared problem, because the bare mass is never equal 
to the physical electron mass. 

2-4- Anomalous Magnetic Moment 

To compute the anomalous moment, we start from the Brodsky-Drell 
formula [28j derived from the spin-flip matrix element of the electromagnetic 



12 



current. In the one-photon truncation, the formula reduces to 



r I 1 j/2+1 r_iv'+ fc ' 

m e J2 / dke? £(f EE ^-^r-^Uk) | (20) 

an J j=0,2 \i'=0k'=j/2 VSfc' 



L l_l\i+k 




Given the explicit expressions © and ([9]) for the wave functions C^, the 
expression for a e simplifies to 



a 

— m, 



j y\l - y)dydk\ ^(-1) W V^ £ (27) 

J W j=0,2 

2 



1 i/ 2 +! f-lV+ fc 



^ ^ + ^ ~ y ^ + ^ ~ m ^(i - y) 

In the sector-dependent case, m 8 is a physical mass but mi is a bare mass. 
The double pole that occurs in the case of the standard parameterization 
is handled in the same way as for the normalization integrals, as discussed 
above. The integrals can be done analytically. 

In the limit where the PV electron mass m\ is infinite, the expression for 
the anomalous moment is 

j y\\-y)dydk\ (28) 

\ 2 

f-ll 



a 


2 2 




m e z 


7T 






X 






\fc=0 



yml + (1 - y)nl + k]_- m\y{\ - y) 



For the sector-dependent parameterization, the product oqZq is just a, and 
the bare mass mo in the denominator is replaced by the physical mass m e . 
To be consistent with [l9j], we eliminate the projection, which does not affect 
the result significantly, and obtain 

a ' = J f 2 * 1 - V)dydkl (29) 



k=0 



ym 2 e + (1 - y)n 2 k + k\ - m 2 e y{l - y) 
13 



as the sector-dependent form of the anomalous moment in the one-photon 
truncation. In the \i\ — > oo, /i — > limit, this becomes exactly the Schwinger 
written as 



result 26 



28 



a 2 

—TO, 
7T 



dydq 2 J{\ - y) 



2i -m 2 



(30) 



The anomalous moment is infrared and ultraviolet safe, and the integral of 
the A; = term yields l/2m 2 . Of course, this limit cannot be taken, because 
the underlying theory would become inconsistent; however, from Fig. HJ we 
can infer that the k — 1 term of fl29|) can be small for sufficiently large 
and that fi can be taken small enough to attain a result near the Schwinger 
result. 

For the standard parameterization, = a, and the normalization is 



1/(1 + e 2 J 2 ), as given in 



The anomalous moment is then 



am 



tt(1 + e 2 J 2 ) 



y 2 (i 




y)dydk]_ 



(31) 



y)vl + k\-mly{l-y) 



This is also plotted in Fig. [TJ For the plotted range of Hi, the normalization 
factor is unimportant. The deviation from the Schwinger result is due to 
the presence of the bare mass in the denominator, and the fact that the 
anomalous-moment integral is sensitive to constituent masses 



18] 



Although the sector-dependent approach, which inserts m e for mo, would 
appear to be a simple remedy, the complications already discussed make this 
choice undesirable. Instead, we can include the physics that does adjust 
the mass by taking into account the self-energy correction that comes from 



29 , which we consider in the next 



the one-electron/two-photon truncation 
section. 

We can also consider the anomalous gravitomagnetic moment, which can 
be computed from a spin-flip matrix element of the energy-momentum ten- 
sor [30]. A result of zero is necessary for any approximation to be deemed 

As shown in 



faithful to the original symmetries of the field theory [31 



30 



whenever the initial and final states can be expressed as Fock-state expan- 
sions with wave functions that depend on boost-invariant, internal momen- 
tum variables, the contributions to the anomalous gravitomagnetic moment 
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are zero, term by term in the Fock expansions. This is precisely the case here, 
for both the standard and sector-dependent parameterizations, with as 
the Fock expansion to be compared with the expansion used for \Vt(P+,P ± )) 
in Eq. (59) of 30], and the net result for this anomalous moment is therefore 
zero, as required. 

3. Self-Energy Contribution 

3.1. Coupled Equations 

We extend the equations (EJ) and (j7|) for the one-electron/one-photon 
sector to include coupling to the one-electron/two-photon sector 

2 i „2 ,,2 i „2 1 

Cg(q) (32) 



M 2 _ m i + q± _ ri + q± 
0--y) y 



6^(-iy^P + [5 s ,±i /2 V£ s (P-q,P) + 6 s>T y 2 U^_ s (P-q,P)] 



_^\a+b^i/ 



abu 



[V:*{P-i-q,P-i)C2i{i,q) 



and add the equation for the new sector 

M 2 _ m l + (qi± + g2±) 2 _ $ + gL _ fif + g% ± 
(l - yi - y 2 ) yi 2/2 

V 1 + fySi" 



C^(q v q 2 ) (33) 
2 E(- X ) a { ^ [VL(P -q.-q^P- g 2 )C^ 



i-2> 



+K,- s (p-g 1 -g 2 ,p-g 2 )cZsk 2 ) 

VLiP-q.-q^P-g^cfs^) 

+ ur a ^ s (P -q_ x -^P- g^S-sig,)} } • 



This last equation can be solved explicitly for the three-body wave func- 
tion C^ ls . Substitution of this into ( 132]) and retention of only self-energy 
contributions for the two-photon intermediate states yields 



M 



2 mt + ql k + q± 



y,q±) 



Iiw{y,q±) 



y,q±) 



(34) 
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with i = 0, 1 and I = 0, 1, 2. Here the coupling to the bare-electron sector is 
written as 



S£ = ^^(-^^i^^.ii^^^-^a+^^C/^-^P-^P)], (35) 

3 

and the new self-energy contributions are written in terms of the integral 

2ttv^ f (-l) a+b+i '£ b e u dq' 

afe^ Wl (1-y-y') y> y 

x PC {P-i-q,L-i) V:, s (P-j-g,P-q) 

+U2(E -q'-<i,p- q')u v aVs {P -q'-q.P-q)]- 



These, combined with the coupled equation (JHJ) for the one-body amplitude, 
constitute the eigenvalue problem when only the self-energy contributions of 
the two-photon states are included. A diagrammatic representation is given 
in Fig. m 

The self-energy integral can be written in simpler form. The change of 
variables x = q' + /q + = y'/ (1 — y) and k± = q' ± + xq± yields 



Jo 



1 dx d 2 k ± mm' - 2^m fl + 



,,, X 7T A m l+ k 'i v 2 b + k i 

a,b u lv l i- x x 

(37) 

which, on use of I , I\, and J defined in Eqs. ffTTl) and ffT2]) . can be written 

as 

?±) = 16tt 2 (-1)*' [mimJoiAt) - 2(m* + m,,)/^) + J(A,)] , (38) 

with 

A^ + (l-y)M 2 -^±^. (39) 

1/ 

For J we still have the identity J(Ai) = AjJo(Aj). The integrals can be 
evaluated analytically, but, for extreme values of the momentum, such as 
momentum fractions on the order of (mo/wii) 2 ~ 10~ 10 , evaluation of the 
analytic form suffers from round-off error due to the finite precision avail- 
able in floating-point calculations. The self-energy is then best computed 
by numerical evaluation of the longitudinal integrals. Notice that I^i in- 
cludes a flavor changing self-energy, where the index i' is different from the 



16 



index i; this is a result of the flavor changing currents used in the interaction 
Lagrangian (1A.1[) . 

The solution for the coupling constant and the one-electron amplitudes 

Zi is 



a± 



Goo + Gn ± a/ (Gqo — Gn) 2 — 4GioGqi 



167r[GooGn — GiqG, 



01 



[Gn — Gqo]/2 =F \/ (Gqo ~~ Gn) 2 — 4GioGqi 



G 



01 



with 



Gu 



-iy 



■[rriimilo - 2 (mi + mi)h + J] 



(40) 
(41) 

(42) 



M 2 - m 2 1 

and the integrals /„ and J defined in flB.6j) . Details of the solution are given 



in 



Appendix B| We again have a as a function of m and the PV masses, and 



we seek the value of mo that yields the physical value of a. The calculation 
requires numerical quadrature of the integrals J„ and J; this can be done 
with the quadrature schemes discussed in [lif and 29 . 



3.2. Anomalous Moment 

The wave functions can be used to compute the anomalous moment. 
Since J does not satisfy any known identity, it cannot be eliminated. There- 
fore, the calculation must be done at finite mi, and the second PV photon 



flavor must be kept [14]]. However, the calculation can still be done [29 



The results are included in Fig. [TJ The value of the PV electron mass 
mi is chosen to be 2 • 10 4 m e , which was found in the case of the one-photon 



truncation to be sufficiently large [l4(. The ratio of PV photon masses //2/V1 



is held fixed at Vz, and /xj is varied. As for the choice of mi, this choice 



was found sufficient in [14|. The results are consistent with perturbative 
QED, showing only variations expected from numerical errors of order 1% in 
calculating the underlying integrals I n and J. 

That the self-energy contribution brings the result so close to the leading 
Schwinger contribution is clear from the following. The dominant contribu- 
tion to the expression (12T)) for the anomalous moment is the j — 0, % = %' = 0, 
and k = k' = contribution to the first term; the other terms are suppressed 
by the large PV masses that appear in the denominators of the wave func- 
tions C^ ks (k). For the dominant term, the denominator, as determined by 
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(IB.4p and (IB.2j) . is essentially the square of 

^4oo + q± 

1-2/ 2/ 
Jooo = 167T 2 [(m2 + A )7 (A ) - 4m / 1 (A )] (44) 



A 00 = m ° + f Io00 + ql + ^±£- m t, (43) 



with 
and 



A = ^ + (l-yK-^l, (45) 

from the expressions in (157)) and (1391) . For the physical photon, /i is zero, 
and the two-body wave function is peaked at q± = and y = 0, so that we 
can approximate Ao as m 2 . We then have 

^ mg + 87ra[(mg + m')/ (m 2 ) -4/720/! (m 2 )] +gj gj 2 

A)o - : 1 rn e (.46) 

1-2/ 2/ 



In our formulation, the perturbative one-loop electron self-energy can be 
read from Eq. (fit))) for i = 0, with Zi = 0, M 2 = itLq + 5m 2 on the left, and 
M 2 = TOq on the right. This yields 

Sm 2 =ml-m 2 e = 2e 2 [m 2 I (ml) - 4m h(ml) + J(mD] . (47) 



From this mass shift, we have, to leading order in a = e 2 /47r, with use of the 
identity J (mo) = mll (ml), 

m 2 = m 2 e - Sixa[(ml + m 2 e )I (m 2 e ) - AmJi(m 2 e )}. (48) 

Therefore, the denominator reduces to the square of 

A 00 = ^±^+ q l-m 2 e + O(a 2 ), (49) 
1-2/ 2/ 

which matches the denominator of the integral ( 130)) that yields the Schwinger 
contribution. Thus, the dominant contribution to the anomalous moment 
with the self-energy included is essentially the same as the integral that 
yields the Schwinger result. 



18 



4. Summary 



In the calculation presented here, we have continued the development 
of nonperturbative Pauli-Villars regularization and light-front Hamiltonian 
techniques as a method for the determination of bound-state wave functions 



in quantum field theories. The new results extend previous work on QED [12 



14 1 to include some of the effects of one-electron/two-photon Fock states; the 



previous work was limited to a one-electron/one-photon truncation. To reach 
agreement with the Schwinger result for the anomalous moment 26( is meant 
as a test of the method but not as the purpose of the method. Instead, the 
method is intended for strongly coupled theories, where perturbation theory 
is not applicable. Applications of the method to a gauge theory are important 
as precursors to applications to QCD. 

We have also shown that use of a sector-dependent parameterization 17, 



18l . Il9| requires great care in the handling of an infrared divergence and an 



uncanceled ultraviolet divergence. The expression for the anomalous moment 
is itself safe from these divergences, and one is led to think that regulators 
can be removed. However, the underlying theory is not safe and quickly 
enters an unphysical regime, if care is not taken. The bare coupling becomes 
imaginary, the Fock-state wave functions become ill-defined, and the Fock- 
sector probabilities fall outside the interval [0, 1]. 

The infrared and ultraviolet divergences are interconnected. The infrared 
divergence is regulated by the introduction of a nonzero mass for the physical 
photon, but the result for the anomalous moment of the electron does not 
agree with experiment unless this mass is quite small. In turn, a small 
physical photon mass requires a large PV photon mass to keep Fock-sector 
probabilities within the physical range of zero to one. Unfortunately, a large 
PV photon mass can lead to great difficulties for a numerical calculation with 
a higher-order Fock truncation, as was seen in 29 . 

The standard parameterization does not have an infrared problem and 
has no particular difficulty with its own uncanceled ultraviolet divergence. 
The Fock-state wave functions are well defined. The one-electron/one-photon 
truncation does yield 14J an anomalous moment that is 17% larger than the 
experimental value; in a sense, it includes too much physics from higher 
orders in a. The discrepancy with experiment, which is unrelated to the 
uncanceled divergence, is immediately corrected by the self-energy contribu- 
tion from the one-electron/two-photon sector, which adjusts the constituent 
electron mass that appears in the anomalous moment integral. The sector- 
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dependent approach, in effect, attempts to incorporate this correction by 
forcing the constituent mass to be equal to the physical mass, which causes 
the infrared divergence. Thus, the standard approach is to be preferred, at 
least for gauge theories. 
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Appendix A. Feynman-gauge QED 

The PV-regulated Feynman-gauge QED Lagrangian is 

2 



i=0 

1 



(A.l) 



i=0 

where 

2 1 

A^E^ 4 ^' ^ = E^' Fi^ = d^A iv - d v A itl . (A.2) 

i=0 i=0 

The subscript i — denotes a physical field and i = 1 or 2 a PV field. Fields 
with odd index i are chosen to be negatively normed. The constants & satisfy 



the following constraints [14 



6 = 1, 

2 



i=0 
2 



i=0 



Mtf/ml) = 0. 



(A.3) 
(A.4) 

(A.5) 



The second constraint guarantees the necessary cancellations for ultraviolet 
regularization; it also implies that A 11 in (1A.2I) is a zero-norm field. The third 
guarantees the correct chiral limit at one loop; for truncations that include 
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higher loops, there are order-a corrections to the constraint. Implementation 
of the gauge condition d^A^ = is discussed in [14 . 

When the PV electron mass is sufficiently large, the third constraint flA.5j) 
can be approximated b}0 



^(-l)V/ 2 ln(/i//m 1 )=0. 



(A.6) 



The solution to the set of constraints, ( IA.3j) . (1A.4j) . and (1A.6j) . assuming 
Ho = 0, is then 



6 = 1 + 6 and 6 



/^ln^/mi) - /x?ln(/xi/mi) 



(A.7) 



Without loss of generality, we require // 2 > so that 6 is positive. 
The dynamical fields are 



A 



1 



V167T 3 



7T° 



7f 



«V(&)e-*« + aL(fc)e*» 



(A.8) 
(A.9) 



with 32[ Xs an eigenspinor of A + = 7°7 + /2. The creation and annihilation 
operators satisfy (anti) commutation relations 



{bUh),blAb!} 

{*.(*), 4y(£} 
[ot/x(fc),aL(^'] 



;-i)'^^5(fc-fc'), 

;-l)^e^(A;-A;')- 



(A.10) 
(A.ll) 
(A.12) 



Here e M = (— 1, 1, 1, 1) is the metric signature for the photon field components 
in Gupta-Bleuler quantization 33|, [34 . 



2 In 



14j , the approximation to this third constraint and the solution for £2 are written 



incorrectly, without the factors of mi and without correct subscripts for the PV photon 
masses, n\ and 
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Without antifermion terms, the Hamiltonian is 

m 2 + p\ 



v 



= £ / 

i,s J 

+£ 



p-T 



l) 4 6t (p)^ s (p) 



(A.13) 



+ ^ I dpdq{bl(p)[baq)V^ 2s (p,q) 

+ h-s(g)utj,-2 S ip,i)] Vti4M-i) + H - c ) ■ 

The vertex functions are [l2| 

eo P± ■ q± ± «P± x <f± + mj-mj + 



\/l6 



7P 



p + q + Vq + — p H 



(A.14) 



-\7l67T 3 



p + q + yjq + — p + 



V* j± (p,q 



Vij±(p,g 



Uf j± (p,q 
Ui i± {p,q 



u ij±(p,g 



For the standard approach, m is the same for all sectors and eo is just e, 
there being no fermion-antifermion loops included in the calculation. For the 
sector-dependent approach, mo and eo depend on the Fock sector where the 
Hamiltonian is applied. The sector-dependent constants could more generally 
be functions of momentum 17], |l8j, but here we follow 19] and use them as 
constants. 



a/167t 3 


P + q + Vq~ 


- -p+ 


e 


p+(q 2 ^fiq 1 ) + 


q + (p 2 ± ip 1 ) 


a/167t 3 


p + q + y/q- 




=Fe 


rrij{p 1 ± ip 2 ) — 


rriilq 1 ± iq 2 ) 


a/167t 3 


p + q + Vq 


f _ p + 


±e 


rrijip 1 ± ip 2 ) — 


rrii(q l ± iq 2 ) 


a/167t 3 


p + q + Vq 




±e 


rriiq + — m,jp + 




Vl67T 3 


p + q + \/q + — p + 




ie 


rriiq + — rrijp + 




-\7l67T 3 


P + q + Vq + ~ P + 
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Appendix B. Solution for the Self-Energy Contribution 



The two-body integral equations (134"]) can be expressed compactly as 
where A^ and Bj are defined by 

Aij = ai±A + ai±A + ±J*_ _ M * m 

' i - v v 2*1 -y 1 



and 



Bj = = -^J^. (B.3) 



We solve this 2x2 system for the two-body wave functions 

A 0j A Xj + B) 



(B.4) 



Without the self-energy contributions, we have Bj = 0, and the wave func- 
tions reduce to the forms given in ([8]) and (Q, with their line of poles whenever 
m < m e , as discussed in Sec. [2j Here, however, the self-energy contributions 
make the denominators more complicated. For values of rriQ that are smaller 
than m e by an amount of order a, there need not be a line of poles. In fact, 
we find that, for the solution with self-energy contributions, there is no pole 
in this two-body Fock sector. 

Substitution of ( IB. 41) into ([5]), and use of the expressions ( 1A.14I) for the 
vertex functions, yields 

[M 2 - ml\zi = 2e 2 ^(-l)^- [m^Jo - 2(m* + m^h + J], (B.5) 
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where 



dydqj ^ . A 0j -A lj -2B j 

167r2 U V Z> y[AojAlj + B]} ' (B - 6) 
dydq]_ \ ^ ( _ j^.'"'- V' ~ m o-4ij - (™o + 



l 6?r 2 , 2/ (i_ ?/ )[A 0j Ai J +^ 2 ] 



X 



16vr 2 

+ gi)A 0i - (m§ + gD^ij - 2(m mi + q\)Bj 



y(l - y) 2 [A^A Xj + B]\ 

When the self-energy contributions are neglected, these return to the previous 
expressions Eq. ( fTTj) and ( |T2l) for 7 , I±, and J in the one-photon truncation. 
What is more, the eigenvalue equation for Zj has nearly the same form as 
the eigenvalue equation (flOj) in the one-photon case. The only difference in 
finding the analytic solution is that Jo an d J are not connected by any known 
identity. 
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Figure 1: The anomalous moment of the electron in units of the Schwinger term (a/2ir) 
plotted versus the PV photon mass, The solid line is the standard-parameterization 
result for the one-electron/one-photon truncation and the dashed lines are the results 
for sector-dependent parameterization at three different values of photon mass (io; both 
use mi = oo. The results for the sector-dependent parameterization are plotted only for 
values of fii for which the probability of the two-particle sector remains between and 
1; for values of < 0.01m e , this requires [i\ > 1000m e , which is beyond the range of 
the plot. For the case of the standard parameterization, /io has its physical value of zero. 
The filled circles arc from a calculation with the standard parameterization that includes 
the self-energy contribution from the one-electron/two-photon sector. It also includes a 
second PV photon flavor, with its mass, /Lt2, set to \/2/ii; the PV electron mass mi is equal 
to 2 • 10 4 m e . The variation is due to errors in numerical quadratures. 
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Figure 2: Diagrammatic representation of the coupled equations §5§ and (|34|) of the text. 
The filled circles and ovals represent wave functions for Fock states; the solid lines represent 
fermions; and the dashed lines represent photons. The crosses on lines represent the light- 
cone kinetic energy contributions, which are summed over all particles in the Fock state. 
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